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Abstract. 

Measurements  of  test-particle  diffusion  in  pure  ion  plasmas  show  2D  enhancements  over  the  3D 
rates,  limited  by  shear  in  the  plasma  rotation  a>£(r).  The  diffusion  is  due  to  “long-range”  ion-ion 
collisions  in  the  quiescent,  steady-state  Mg+  plasma.  For  short  plasma  length  Lp  and  low  shear 
S = r dotE/dr , thermal  ions  bounce  axially  many  times  before  shear  separates  them  in  0,  so  the 
ions  move  in  (r,  0)  as  bounce  averaged  “rods”  of  charge  (i.e.  2D  point  vortices).  Experimentally,  we 
vary  the  number  of  bounces  over  the  range  0.2  <Nb<  10,000.  For  long  plasmas  with  Nb  < 1,  we 
observe  diffusion  in  quantitative  agreement  with  the  3D  theory  of  long-range  E x B drift  collisions. 
For  shorter  plasmas  or  lower  shear,  with  Nf,  > 1,  we  measure  diffusion  rates  enhanced  by  up  to 
100 x.  For  exceedingly  small  shear,  i.e.  Nb  > 1000,  we  observe  diffusion  rates  consistent  with  the 
Tayior-McNamara  estimates  for  a shear-free  thermal  plasma.  Overall,  the  data  shows  fair  agreement 
with  Dubin’s  new  theory  of  2D  diffusion  in  shear,  which  predicts  an  enhancement  of  D2D/D3D  « Nb 
up  to  the  Tayior-McNamara  limit. 


In  turbulent  fusion  plasmas,  experiments  [1]  have  observed  reduced  transport  in  the 
presence  of  sheared  flow,  but  the  transport  is  difficult  to  calculate  theoretically  [2].  Here 
we  report  measurements  of  test  particle  diffusion  in  a quiescent  magnesium  ion  plasma 
for  direct  comparison  to  a new  theory  of  coliisional  diffusion  in  a 2D  point-vortex  gas 
with  background  shear.  The  ions  are  kept  in  steady  state  with  a rotating  wall  [3]  drive 
superimposed  on  the  confining  potential  at  one  end.  We  observe  that  for  short  plasmas 
(2D)  the  diffusion  rate  is  up  to  100  times  larger  than  for  long  plasmas  (3D).  This  can 
be  easily  understood  if  we  introduce  Nb,  the  number  of  bounces  that  a particle  makes, 
in  the  rotating  frame,  before  being  sheared  apart  (in  the  0 direction)  from  a neighboring 
particle. 

We  find  that  Nb  controls  the  2D  transport  enhancement  over  the  3D  case.  Stated 
differently,  if  Nb  is  large,  neighboring  particles  interact  for  longer  periods  of  time, 
leading  to  more  transport.  For  large  Nb  the  ions  can  be  described  as  bounce  averaged 
“rods”  of  charge  (which  theorists  refer  to  as  “2D  point  vortices”).  In  the  absence  of 
shear,  the  transport  is  limited  by  the  Tayior-McNamara  paradigm  [4],  where  thermal 
fluctuations  excite  large  scale,  but  short  lived,  convective  cells.  In  the  presence  of  shear, 
new  quantitative  theoretical  work  [5]  predicts  that  the  diffusion  coefficient  is  inversely 
proportional  to  the  shear. 

In  our  experiment,  the  transport  is  controlled  by  the  bounce  averaged  shear,  i.e.  the 
average  shear  that  a particle  experiences  averaged  over  one  bounce  along  a magnetic 
field  line  of  the  trap.  To  experimentally  determine  the  E x B shear,  we  use  a laser  induced 
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fluorescence  (LIF)  diagnostic,  giving  the  density  n(r),  the  temperature  T(r)  and  the  fluid 
rotation  velocity  v10t(r)  in  the  G direction  using  a perpendicular  laser  beam  as  shown  in 
Fig.  1.  All  the  LIF  measurements  are  made  in  the  central  plane  (z  = 0)  of  the  plasma.  For 
the  present  work,  we  have  about  108  magnesium  ions,  with  a plasma  length  L,,  ~ 1 cm, 
and  a plasma  radius  rp  ~ 2 cm;  the  temperature  is  controlled  from  0.1  eV  to  4 eV;  and 
the  magnetic  field  is  B = 3 Tesla.  More  details  about  the  apparatus  and  diagnostic  can 
be  found  in  Ref.  [6]. 

The  density,  temperature  and  rotation  profiles  are  shown  in  Fig.  2 for  a plasma  with 
small  shear  (a)  and  for  a plasma  with  large  shear  (b).  We  fit  the  rotation  profile  vtl)l  to 
a smooth  third  order  polynomial.  We  numerically  solve  Poisson’s  equation  using  the 
smooth  fit  to  V[0(  and  T(r)  with  the  wall  potentials  as  boundary  conditions;  and  obtain 
n(r,z)  and  <t>(r,z).  From  this  solution,  one  obtains  the  shear  rate  S(r,z)  = rdwi;(r,z)/dr 
and  the  bounce  averaged  dimensionless  shear  rate  [5] 


(s{r))z 


S{r,z)  \ 
2nn(r,z)  ce/B  / z 


The  dimensionless  shear  rate  s(r)  represents,  for  the  case  of  a “top  hat”  profile,  the  shear 
rate  divided  by  the  local  rotation  rate.  Figure  2 also  shows  (s)~\  the  squares  represent 
the  bounce  average  shear,  the  solid  line  a simple  estimate  of  s(r,z  = 0)  directly  from 
v £■(/•)  = V(0((r)  — V(jia(r)  where  the  diamagnetic  drift  is  calculated  from  the  measured 
/i(r)  and  T(r).  In  the  case  of  Fig.  2(a)  which  has  a very  small  shear  for  small  radius, 
(s(r))z  ~ .r(/-,z  = 0)  for  r < 1 cm. 

The  electronic  spin  orientation  of  the  ground  state  of  Mg+  is  used  to  “tag”  the  test 
particles.  The  flux  of  test  particles  T,  is  obtained  from  the  measured  test  particle  density 
n,  as 


T ,(r,t)  = --J  dx x~n,(x,l)  + J dxx 


2n,(x,t)  — n(x) 

Tv  (a‘) 
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FIGURE  2.  Profiles  of  density  n,  temperature  T,  velocity  v,  and  dimensionless  shear  s.  (a)  Small  shear 
case;  (b)  large  shear  case. 


the  second  term  corrects  for  the  weak  spontaneous  spin  flip  at  a measured  rate  t5-1  . We 
verify  that  the  flux  of  test  particles  is  proportional  to  the  gradient  of  the  concentration  of 
test  particles,  i.e. 


r,(r,0  = -£>(r)n(r) 


d_ 

dr 


n(r)  ) 


and  obtain  the  local  diffusion  coefficient  D(r). 

Figure  3 shows  preliminary  measurement  of  the  diffusion  coefficient  normalized  to 
the  plasma  length  Lp  and  magnetic  field  B as  a function  of  the  dimensionless  bounce- 
averaged  shear  (s)z.  For  (s)z  >0.1  the  measured  diffusion  coefficient  agrees  quantita- 
tively wilth  the  3D  diffusion  coefficient  previously  measured  in  long  plasmas  [7];  here, 
D3D  is  calculated  with  values  of  n and  T appropriate  to  the  high-shear  data  points.  As 
the  shear  decreases,  the  diffusion  coefficient  increases  roughly  proportional  to  s.  For 
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FIGURE  3.  Diffusion  coefficient  normalized  to  plasma  length  L,,  and  magnetic  field  II  versus  the 
dimensionless  bounce  average  shear. 

(s)z  < 5 x 10~3  the  diffusion  coefficient  approaches  a limit  within  a factor  of  two  of 
the  Taylor-McNamara  shear-free  theory.  Here,  D™  is  calculated  from  the  number  of 
particles  N in  the  low-shear  plasmas. 

Recent  theoretical  work  has  investigated  diffusion  processes  in  the  presence  of  shear 
for  a collection  of  2D  point  vortices  [5].  Overall,  the  data  shows  fair  agreement  with 
the  present  theory,  but  proper  consideration  of  the  end  shape  and  Debye  shielding  of  the 
fluctuations  will  have  to  be  included.  To  summarize  the  theories,  we  have  the  3D  long 
range  collisional  diffusion  [8 j: 

D30  = 2a.y/nn'jb2  r2  ln(v/Avmjn)ln(Xo/rc.)  °c  B~2iiT] /2, 

where  a = 3 represents  the  multiple  effective  collisions  due  to  “velocity  caging.”  This 
is  in  quantitative  agreement  with  the  experiments  [7]  in  long  plasmas  (i.e.  Nj,  < 1).  Here 
b = e2 /T  is  the  classical  distance  of  closest  approach,  and  the  minimum  relative  velocity 
for  these  collisions  is  Avmjn  ss  (nv3b2  x/rvAy)'/3. 

The  2D  shear-limited  diffusion  [5]  is 


where  d = (4D2D/|5|)'/2  is  the  (diffusion-limited)  minimum  impact  parameter  for  these 
collisions.  In  the  limit  of  zero  shear,  Taylor  and  McNamara  calculated  [4]  a diffusion 
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FIGURE  4.  The  diffusion  coefficient  normalized  by  the  3D  long  range  E x B drift  collision  case 
increases  with  Afc. 


(neglecting  considerations  of  Debye  shielding)  of 


D™ 


ce  \/W 
27t‘/2  bTp 


These  theories  scale  differently  with  all  plasma  parameters.  However  the  enhancement 
ratio  D2D/D3D  depends  only  on  the  number  of  bounces  a particle  makes  in  the  rotating 
frame  before  being  sheared  apart,  given  by 


Nb  = 


Jb_ 

m'E 


v 

2Lpru)'E 


Figure  4 shows  the  measured  diffusion  coefficient  normalized  to  D3D.  The  solid  line  is 
the  theory  prediction  at  B = 3 T for  D2D/D3D  « Nj,.  One  can  see  that  the  diffusion 
increases  as  particles  interact  for  longer  periods,  i.e.  for  large  Nf,.  For  Nb  <2  the 
diffusion  is  correctly  described  by  3D  long-range  collisions;  in  contrast  for  Nb  > 5000 
the  diffusion  approaches  the  Taylor-McNamara  limit.  For  \ <Nb  < 100,  the  predicted 
D2D/D3D  is  in  fair  agreement  with  the  experiments,  being  generally  about  3x  larger 
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than  observed.  Further  theory  and  experiments  may  elucidate  the  effects  of  the  end 
confinement  fields  and  the  effects  of  Debye  shielding. 
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